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1.  Projectile  Dynamic  Model 


This  report  extends  the  work  of  Murphy  and  Bradley1  and  Fansler  and  Schmidt2  describing  jump 
caused  by  slight  configuration  asymmetries.  Some  of  their  results  are  briefly  repeated  here,  for 
completeness,  followed  by  a  presentation  of  a  closed-form  analytic  solution.  The  discussion  then 
continues  with  a  simple  rational  approximation,  which  is  then  compared  to  the  analytic  solution. 


Flight  mechanics  of  most  projectile  configurations  can  be  captured  using  a  rigid  body  six  degrees 
of  freedom  dynamic  model.3  The  degrees  of  freedom  are  three  position  components  of  the 
projectile  mass  center  and  three  Euler  orientation  angles  of  the  body.  Figures  1  and  2  show  two 
helpful  schematics  so  that  the  degrees  of  freedom  are  seen  to  be  related  according  to  the 
following  equations  of  motion,4 
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1  Murphy,  C.  H.;  Bradley,  J.  W.  Jump  Due  to  Aerodynamic  Asymmetry  of  a  Missile  With  Varying  Roll  Rate ;  BRL-R-1077; 
U.S.  Army  Ballistic  Research  Laboratory:  Aberdeen  Proving  Ground,  MD,  1959. 

2 

“  Fansler,  K.  S.;  Schmidt,  E.  M.  Trajectory  Perturbations  of  Asymmetric  Fin-Stabilized  Projectiles  Caused  by  Muzzle  Blast. 
Journal  of  Spacecraft  and  Rockets  1978, 15  (1),  62-64. 

3 

McCoy,  R.  L.  Modern  Exterior  Ballistics:  the  Launch  and  Flight  Dynamics  of  Symmetric  Projectiles',  Schiffer  Publishing 
Ltd.:  Atglen,  PA,  1999. 

4  Murphy,  C.  H.  Free  Flight  Motion  of  Symmetric  Missiles;  BRL-TR-1216;  U.S.  Army  Ballistic  Research  Laboratory: 
Aberdeen  Proving  Ground,  MD,  1963. 
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Figure  1.  Projectile  position  coordinates  definition. 


Forces  in  the  body  frame  that  appear  in  equation  3  contain  contributions  from  weight  (W)  and  air 
loads  (A), 
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The  weight  force  resolved  into  projectile  body  coordinates  is  given  by  equation  6, 
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The  air  loads  are  split  into  two  components:  the  standard  aerodynamic  forces  and  the  Magnus 
forces, 
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Figure  2.  Projectile  orientation  definition. 


Equation  8  gives  the  standard  air  loads  acting  at  the  aerodynamic  center  of  pressure, 

(v2+w2) 

I  Wo  _r'-'X2 

N 

>  =  -q*  < 


Y 


AS 


-'AS 


V2 


c  — 

NA  V 


c  - 

NA  V 


where 


and 


qa 


-p(u2  +  v2  +  w2)7rD2 
8 


V  =  -\/u2  +  v2  +  w2  . 


The  Magnus  aerodynamic  force  acts  at  the  Magnus  center  of  pressure, 
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Moments  about  the  projectile  mass  center  are  due  to  aerodynamic  forces,  and  moments  (A)  are 
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The  aerodynamic  moments  caused  by  standard  and  Magnus  air  loads  are  computed  with  a  cross 
product  between  the  distance  vector  from  the  mass  center  to  the  force  application  point  and  the 
force  itself.  An  unsteady  aerodynamic  damping  moment  is  also  present,  which  provides  a 
damping  source  for  angular  motion, 
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All  aerodynamic  coefficients  and  the  center  of  pressures  are  a  function  of  the  Mach  number  of 
the  projectile  mass  center.  The  dynamic  model  previously  described  is  nonlinear  due  to  both 
three-dimensional  rotational  kinematics  expressions  and  the  presence  of  complex  aerodynamic 
forces. 


2.  Reduction  to  Linear  Theory 


Useful  performance  data  regarding  trajectory  prediction  and  the  stability  of  projectiles  forced 
early  ballisticians  to  investigate  mathematical  simplifications  to  the  equations  of  motion.  Over 
time,  a  set  of  simplified  and  solvable,  yet  accurate,  linear  differential  equations  emerged,  which 
today  is  commonly  tenned  “projectile  linear  theory.” 

The  governing  equations  previously  developed  are  expressed  in  the  body  reference  frame.  In 
linear  theory,  the  lateral,  translational,  and  rotational  velocity  components  are  transformed  to  a 
nonrolling  reference  frame.  The  nonrolling  frame,  or  so-called  fixed  plane  frame,  proceeds  with 
only  precession  and  nutation  rotations  from  an  inertial  reference  frame.  Components  of  linear 
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and  angular  body  velocities  in  the  fixed  plane  frame  can  be  computed  from  the  body  frame 
components  of  the  same  vector  through  a  single-axis  rotational  transformation.  For  example,  the 
body  frame  components  of  the  projectile  mass  center  velocity  are  transformed  to  the  fixed  plane 
frame  by 
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Note  that  the  :  superscript  indicates  the  vector  components  relative  to  the  fixed  plane  reference 
frame.  Projectile  linear  theory  makes  a  change  of  variables  from  station  line  velocity 
component,  u,  to  total  velocity,  V,  as  described  in  the  next  two  equations: 

V  =  Vu2  +  v2  +  w2  =  Vu2  +  v2  +  w2  ,  (15) 
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A  further  change  of  variable  from  time,  t,  to  dimensionless  arc  length,  s,  is  also  preferred  and 
following  Murphy4  gives  the  dimensionless  arc  length, 
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(17) 


Equations  18  and  19  relate  time  and  arc  length  derivatives  of  a  given  quantity^  .  Dotted  tenns 
refer  to  time  derivatives,  and  primed  terms  denote  dimensionless  arc  length  derivatives, 
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Linear  theory  makes  several  assumptions  regarding  the  relative  size  of  different  quantities  to 
further  simplify  the  analysis.  Euler  angles  are  small  so  sin(9)  «  9  ,  cos(9)  «  1 ,  sin(\j/)  «  \| / ,  and 

cos(v| /)  «  1 ,  and  the  aerodynamic  angle  of  attack  is  small  so  that  a  =  w/V  and  (3  =  v/V .  The 
projectile  is  mass-balanced  such  that  IXY  =  Ixz  =  IYZ  =  0  and  Izz  =  IYY  =>  IYY  =  Izz  =  IY  . 
Quantities  V  and  <f>  are  large  compared  to  9,  v| /,  q,  r,  v  ,  and  w  ,  such  that  products  of  small 
quantities  and  their  derivatives  are  negligible.  Application  of  these  assumptions  results  in 
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Aside  from  the  fact  that  V  appears  in  some  of  the  previously  mentioned  coefficients,  the 
dynamics  are  now  expressed  with  linear  ordinary  differential  equations. 


3.  Linear  Theory  Solution 


Linear  theory  offers  physical  insight  into  the  flight  dynamics  because  closed-form  solutions  can 
be  readily  obtained.4  Because  V  changes  slowly  with  respect  to  the  other  variables,  it  is  thus 
considered  constant,  V  «  V0 ,  when  it  appears  as  a  coefficient  in  all  dynamic  equations  except  its 

own.  Moreover,  pitch  attitude  of  the  projectile  is  regarded  as  constant  in  the  velocity  equation, 
thus  decoupling  the  velocity  equation.  The  epicyclic  motion,  equation  28,  together  with  the  roll 
dynamics,  equation  27,  is  uncoupled  and  forms  a  linear  system  of  equations.  In  projectile  linear 
theory,  the  Magnus  force  in  equations  24  and  25  is  typically  regarded  as  small  so  that  in  further 
manipulation  of  the  equations,  all  Magnus  forces  will  be  dropped.  However,  it  is  important  to 
retain  Magnus  moments  due  to  the  fact  that  a  cross  product  between  Magnus  force  and  its 
respective  moment  arm  is  not  necessarily  small. 


The  solution  to  the  differential  equation  26,  for  the  forward  velocity,  is 


(  — 2ayS 

\ 

V0e  v  +-*-1 

e 

-1 

J  av 

V 

) 

(30) 


When  90  =  0 ,  the  velocity  solution  reduces  to  the  familiar  exponential  decay  form.4  The  roll 

dynamic  equation  is  a  nonhomogeneous  linear  differential  equation  with  the  following  solution 
for  90  =  0 : 
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Noting  that  p(s)  =  ^'V/D  implies 
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Neglect  the  product  of  damping  and  the  product  AE  since  the  density  ratio  is  assumed  small. 
Defining  £  =  (v  +  iw)/V0  enables  equation  28  to  be  reduced  to  a  single  differential  with  the 

following  form,4 

4"-(E- A)£'-(AE  +  C)4  =  -(MM  +  EFF)e'^+^  -iEG/V0 ,  (33) 

by  assuming  p  is  small  and  ignored. 
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Rather  than  solve  for  the  lateral  translation  and  rotational  velocity  components,  via  equation  33, 
a  more  direct  way  to  obtain  the  effects  of  asymmetry  is  to  solve  the  swerve  differential 
equations.  The  lateral  translation  and  rotational  velocity  components  are  contained  in  the 
attitude  differential  equations,  and  the  attitudes  are  contained  within  the  swerve  differential 
equations. 


4.  Swerve 


Swerving  motion  along  the  earth-fixed  J,  and  K,  axes  results  from  a  combination  of  the  nonnal 

aerodynamic  forces,  as  the  projectile  pitches  and  yaws,  plus  the  forces  and  moments  due  to  the 
configuration  asymmetry.  Differentiating  equations  22  and  23,  with  respect  to  nondimensional 
arc  length  and  using  the  definition  of  £,  with  equation  33,  leads  to  the  following  expression, 

(AE  +  C)  +  iz^  =  _ (AMM  _  CFF)  i  +  (E  -  A)A5'  -  A^" .  (34) 

D  V0 

For  a  stable  projectile,  the  swerve  caused  by  epicyclical  vibration  decays  as  the  projectile 
progresses  downrange  and  does  not  affect  the  long-tenn  lateral  motion.  However,  the 
assumption  that  the  projectile  is  configurationally  asymmetric  causes  an  integrated  effect  that 
contributes  to  the  long-term  lateral  motion  of  the  projectile.  Linear  theory  shows  this  center  of 
mass  motion  contains  terms  that  are  bounded  with  arc  length  s  plus  terms  that  are  linear  with  s 
and  with  the  inclusion  of  gravity  the  solution  of  equation  34  will  have  even  higher  order 
diverging  terms.  These  higher  order  terms  are  typically  denoted  as  gravity  drop.  The  linear 
tenns  are  called  jump  terms,  which  are  caused  by  initial  conditions  at  the  gun  muzzle,  forces 
caused  by  asymmetry,  and  aerodynamic  characteristics.  Ignoring  gravity  and  evaluating  the 
following  limits  formally  defines  aerodynamic  jump 

lim^  =  TJ5  and  lim^  =  TK .  (35) 

S-HO  D  S  S-»oo  D  s 

The  total  aerodynamic  jump  vector  T  is  expressed  as  the  sum  of  two  vectors.  The  first  vector 
represents  the  muzzle  conditions  and  the  second  results  from  asymmetry  subjected  to  a  varying 
roll  rate: 
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The  quantity  II  is  the  contribution  to  the  jump  vector  attributed  to  the  assumed  asymmetry  of 
the  projectile.  The  appendix  shows  that 


n  = 


AMM-CFF 
AE  +  C 


i  +  I^Tz  ■F.(1;2-iz'“irz) 


(38) 


where  |  F  is  the  confluent  hypergeometric5  function.  Apparently,  neither  Murphy  and  Bradley1 

nor  Fansler  and  Schmidt”  were  aware  of  equation  38,  which  makes  their  limiting  and  asymptotic 
analysis,  T  =  1  and  |r  - 1|«  1  ,  unnecessary  for  calculations  of  jump  due  to  asymmetric 

configurations. 


For  constant  rolling  motion  T  =  0  and  for  comparison  purposes,  the  following  definition 
introduced  by  Murphy  and  Bradley1  will  be  used  here,  FI  =  -ill  Z ,  so  that  O  =  -iOZ  of 
equation  (A-6)  of  the  appendix  now  becomes 

6  =  l-irzje"z(ry_ln(1_y))dv 
0 

=  l-iTZ(l-iZ)  1  jFj  (1;2  -  Zi, - iTZ).  (39) 


Polar  plots  of  results  from  equation  39  composed  to  Murphy  and  Bradley1  are  given  in  figure  3, 


for 


<D 


,  and  figure  4  gives  the  argument  O  from  equation  39. 


Magnitude  of  <j> 


Figure  3.  Magnitude  of  O  as  a  function  of  Z . 


5  Abramowitz,  M.;  Stegun,  I.  A.,  Eds.  Handbook  of  Mathematical  Functions,  Series  55;  National  Bureau  of  Standards 
Applied  Mathematics:  Dover  Publishing  Co.,  1967. 
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Figure  4.  Argument  of  O  as  a  function  of  Z  . 


Successive  integration  by  parts  of  the  last  equation  shows 
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Hence,  in  the  limit  of  large, 
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Furthermore,  equation  39  can  be  expanded  using  a  Kurnmer  series  giving  the  following 
expression: 
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(42) 


The  limit  of  equation  42,  Z  — »  oo ,  is 

6  =  i+r  +  r2 +r3 +r4---  =  forr<i. 


(43) 
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This  last  result,  in  light  of  the  general  limiting  case  of  equation  41,  suggests  the  Kummer  series, 
equation  42,  when  transformed  to  a  continued  fraction,  may  produce  an  accurate  approximation 
to  O  for  all  values  of  T  .  To  continue  the  investigation,  the  assumption  was  made  that  a 
reasonable  approximation  has  the  following  representation: 

6  =  aOT/(l  +  air/(l  +  a2T/(l  +  a3r /• ,  -(l  +  a7T ))))))) .  (44) 

After  finding  the  coefficients  aO,  al,  a2---a7  results  in  the  rational  approximation 


®  = 


iZ5  (Z  +  2i)(Z  +  3i)(Z  +  4i)(Z  +  5i)2  (Z  +  6i)(Z  +  7i)(z(Z  +  5i)  -  10)r4 

-  iZ3  (Z  +  2i)(Z  +  3i)(Z  +  4i)2  (Z  +  5i)2  (Z  +  6i)2  (Z  +  7i)(z(3Z  +  5i)  -  8)r3 
+  3iZ3  (Z  +  2i)(Z  +  3i)2  (Z  +  4i)2  (Z  +  5i)2  (Z  +  6i)2  (Z  +  7i)2  T2 

-  iZ(Z  +  2i)2  (Z  +  3i)2  (Z  +  4i)2  (Z  +  5i)2  (Z  +  6i)2  (Z  +  7i)2  (Z  +  8i) T 
i(Z  +  i)(Z  +  2i)2  (Z  +  3i)2  (Z  +  4i)2  (Z  +  5i)(Z  +  6i)(Z  +  7i)Z 

Z3r4  -4Z(Z  +  5i)r3  +6Z(Z  +  5i)(Z  +  6i)r2  -4(z  +  5i)(z  +  6i)(z  +  7i)r+  ^ 
_(Z  +  13i)(Z(Z  +  13i)-82) 


•  (45) 


Taking  the  limit  of  this  expression,  for  Z  — >  oo ,  yields  ®  — > 


1 

1-T 


which  indicates  equation  45  is 


a  reasonable  approximation  when  T  ^  1 .  Figure  5  shows  some  comparisons  between  the  exact 
solution,  equation  39,  and  its  rational  approximation  equation  45  for  various  values  of  T .  It  is 
noteworthy  that  even  for  the  case  where  T  =  1 ,  the  comparison  shows  good  agreement  for 
Z  <  20 . 


Figure  5.  A  comparison  of  the  exact  vs.  a  rational  approximation  of  ®  . 
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5.  Conclusions 


Previous  efforts  describing  jump  due  to  asymmetry  did  not  develop  the  closed-form  solution 
presented  here.  All  of  the  analysis  is  based  on  projectile  linear  theory,  which  leads  to  an 
expression  based  on  the  confluent  hypergeometric  function  iF i .  This  solution  is  well 
approximated,  for  the  arguments  used  here,  with  a  simple  rational  expression  obtained  from  a 
continued  fraction  expansion  of  the  closed-form  solution.  This  is  a  further  extension  of  the 
Murphy  and  Bradley  and  Fansler  and  Schmidt  results  (see  figure  3)  that  will  prove  useful  for 
analysis  and  design  purposes. 
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Appendix.  Integral 


Let  f  =  1  -  (/)'(]  / (f>[  ,  Z  =  (/)[  / K p  and  x  =  sKp ,  then  equation  32  will  be  written  as 


^(t)  =  Z  x  +  T  e  T-1  . 


(A-l) 


Then 


©  =  hm^-jje^^dc  dx 

s  0  0 


limj^oo  ~~J (s  ~  x)  e' 


Z  T  +  r(e_,-l)] 


(A-2) 


Ignoring  the  limiting  process  for  the  moment  allows  the  last  equal  sign  to  be  written  as 

s®  =  +r  £  IZLLj (S  -  i)e|lz -”>'dT . 


(A-3) 


Integrating  the  last  expression  and  taking  the  limit  s— >co  causes  equation  A-3  to  become 


0  =  e  izr  - 


i  |  y  (lZr)n  ' 

Z  n4l(n-iz)nlj’ 


(A-4) 


where  after  writing  the  summation  as  two  sums  over  even  and  odd  values  of  n  respectively 
becomes 


0  =  e'izr  + 


(-ir(zr)^  , 

■-1  (2n)!  I 


(A-5) 


„  (-i)n(zr)2n+1 j.  .... 

+i  y - — aa — 4 —  f  t  "  ]  dt 

“0  (2/7  +  1)!  J 


Using  Taylor  Series  expansions  followed  with  partial  integration  given  by  A  and  S  , 
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=  -+  r  ^(va-iz,  -  izr) 

Z  1  -  iZ  1  u  ’ 


(A-6) 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


Integration  constant  a 


pClddD  S 


Integration  constant  av 


pCxoDS 


Integration  constant  bp  = 


pCLp  D3  S 


Integration  constant  bv  = 


g90D 


Projectile  aerodynamic  coefficients. 
Projectile  characteristic  length  (diameter). 
Force  component  caused  by  asymmetry. 
Gravitational  constant. 

Scaled  gravitational  constant  G  =  gD/V0 


Mass  moments  of  inertia. 


Constant  defined  as  K  =  - 


pDS(2CxoII  +  mCLPD!) 


Applied  moments  about  projectile  mass  center  expressed  in  the  no-roll  frame. 


Moment  component  caused  by  asymmetry. 


Projectile  mass. 


Angular  velocity  components  vector  of  projectile  in  the  no-roll  frame. 


Dummy  integration  variable. 
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s 

u 

V 

w 

V0 


Z 

a 

P 

r, 

rK 

r 

n 

4> 

e 

¥ 

0  0 
£ 

P 


Surface  area  S  =  tiD  /4. 


Mass  center  velocity  components  in  the  body  reference  frame. 


Forward  velocity  of  projectile. 


Mass  center  velocity  components  in  the  no-roll  reference  frame. 


Position  vector  of  body  center  of  mass  in  an  inertial  reference  frame. 


Substitution  variable  defined  as  Z  =  ^'/Kp  . 

Longitudinal  aerodynamic  angle  of  attack. 

Lateral  aerodynamic  angle  of  attack. 

' 1  Components  of  aerodynamic  j  ump . 

K; 

Substitution  variable  defined  as  T  =  I  -  (f)'t]  / (f)[  . 
Complex  aerodynamic  jump  caused  by  asymmetry. 

Euler  roll,  pitch,  and  yaw  angles  of  the  projectile. 


Euler  roll  rates 


<K 

€ 


Po  d/v, 

aPD/V 


Euler  roll  angle  of  configuration  asymmetry. 


Air  density. 
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